NMR linewidth and Skyrmion localization in quantum Hall ferromagnets 
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The non-monotonic behavior of the NMR signal linewidth in the 2D quantum Hall system is 
explained in terms of the interplay between skyrmions localization, due to the influence of disorder, 
and the non-trivial temperature dependent skyrmion dynamics. 
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I. INTRODUCTION 



The ground state of the two-dimensional electron gas 
(2DEG) at filling factor one [y = 1) is a spin-polarized 
stated in which all electrons completely fill the lowest 
Landau level with spin up polarization (quantum Hall 
ferromagnet) . The low lying charged excitation of the 
2DEG is called skyrmion, or charged spin texture, and 
carries an unusual spin distribution^. Pictorially, this dis- 
tribution can be viewed as a configuration in which the 
spin field points down at a given position and smoothly 
rotates as one moves radially outwards from that point, 
until all the spins are polarized as in the ground state. 
These non-trivial objects are topologically stable, with its 
size -the number of reversed spins- fixed by the compe- 
tition between the Coulomb and Zeeman interactions 3 ;. 
The existence of this many body electronic state was 
confirmed by Barrett and coworkers^ who performed op- 
tically pumped nuclear magnetic resonance (OPNMR) 
measurements of the 71 Ga nuclei located in an n-doped 
GaAs multiple quantum well. As the nuclei of the quan- 
tum wells are coupled to the the spins of the 2DEG via 
the isotropic Fermi contact interaction, the Knight shift 
{Kg) from the NMR signal gives information about the 
spin polarization of the 2DEG. It was reported that K s 
symmetrically drops around v = 1, showing that the 
charged excitation of the quantum Hall ferromagnet in- 
volves a large number of spin-flips^ in contrast with the 
expected behavior for non-interacting electrons. Since 
then, heat capacity, magnetoabsorption and magneto- 
transport measurement s^"^ have been performed in or- 
der to understand the main features of this novel state, 
however, valuable information can be still extracted from 
the NMR experiments. 

Recently, new OPNMR measurements^ of the 71 Ga 
nuclei in n-doped GaAs/Alo.iGao.gAs multiple quantum 
wells were performed at temperatures (T ~ 0.3 K) lower 
than the ones considered in Ref. 0- In Ref. 0, Khan- 
delwal et al. observed that the Knight shift data around 
v = 1 consist of a "tilted plateau", in contrast to the 
symmetrical fall at higher temperatures. As NMR is a 
local probe, it was argued that the existence of the tilted 
plateau is related to the localization of the skyrmions. 
In addition, the full width at half maximum (FWHM) 



of the OPNMR signal has a non-monotonic temperature 
dependence at v w 1, increasing when the temperature is 
lowered up to a critical value and then, falling down. It 
is well knownS that the linewidth of the NMR measure- 
ments is affected by the dynamics of the nuclei {motional 
narrowing effect). However, as the 71 Ga nuclei are spa- 
tially fixed in the GaAs quantum wells, any observed dy- 
namical effect must be related to the dynamics of 2DEG. 
In particular, if v w 1, those effects should be related 
to the skyrmion dynamics. In this scenario, the behav- 
ior of the FWHM can be understood as an evolution of 
the skyrmion dynamics from the motional-narrowed to 
the frozen regimes, as the temperature decreases^. This 
non-monotonic behavior of the linewidth was previously 
observed by Kuzma et aliM when the filling factor of the 
2DEG was slightly changed around v — 1/3. Later on, 
in order to investigate the evolution of the spectral line 
shape as function of T, Sinova and coworkers^ simulated 
the NMR signal of 71 Ga nuclei when i/wl, As at low T 
it is expected that the short length scale correlations con- 
tain the ground state features, the 2DEG in Ref. Hll was 
described by an square skyrmion latticed that randomly 
jumps -as a whole- over an average distance. Despite the 
fact that the formation of a pinned Wigner crystal has 
never been observed in the range of temperatures where 
the experiments&ifi have been performed, the approach 
used in Ref. ^ reproduced the qualitative behavior of 
the FWHM near the completely frozen situation. How- 
ever, as far as we know, no explicit correspondence has 
been established between the temperature dependence of 
the skyrmion dynamics, for the whole range of experi- 
mentally available temperatures, and the line shape of 
the NMR experiments. This is precisely the main goal of 
this paper. 

Here, instead of trying a fully quantum mechanical ap- 
proach for the Fermi contact interaction between 71 Ga 
nuclei and the 2DEG for v w 1 we will use a semiclassical 
theory, in which the skyrmion center of mass is treated 
as a true dynamical variable 13 to which a temperature 
dependent diffusion constant is associated^. This way 
of facing the problem will allow us to study the whole 
range of temperature, namely from the complete delocal- 
ized situation to the frozen limit. 

We organize this paper as follows. In section II the 
model used to describe the skyrmion dynamics and the 
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localization by disorder are presented and the dependence 
of the linewidth with the temperature is calculated. Sec- 
tion III is devoted to analyze the change in the energy 
of the skyrmions when impurities are taken into account. 
Finally in Sec. IV we present the final results and our 
conclusions. 



II. THE LINEWIDTH TEMPERATURE 
DEPENDENCE 

To start with we recall that the FWHM in the mag- 
netic resonance experiment is proportional to the in- 
duced phase shift experienced by the nuclei every time 
they interact with a fluctuating field. More precisely, the 
FWHM is equal to T 2 _1 , where T2 is the transversal relax- 
ation timA For the Ga nuclei in the quantum wells, the 
fluctuating field is given by the deviation from the com- 
pletely polarized state provided by moving skyrmions. 
Therefore, the spins of the 71 Ga nuclei will precess with 
an extra phase over its normal precession (in the presence 
of the completely polarized state) every time a skyrmion 
passes near by. This change in phase can be calculated 
as S(j) = ± 7„ n t— , where j n is the giromagnetic ratio for 
the 71 Ga nuclei, n is the "fluctuating" field associated to 
the skyrmion magnetization and r is the interaction time. 
The mean square dephasing (A0 2 ) after n r-intervals of 
time will be n (j n Tlt) 2 . Here n = 2DtNf s /S is the num- 
ber of skyrmions passing by the nucleus in a time t, with 
Nf s being the number of free skyrmions, D the diffu- 
sion constant and S the quantum well surface area. As 
a result, 
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Now, if we define T2 as the time for a group of spins in 
phase at t = get about one radian out of step& we find 
that 
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On the other hand, the time interval r can be estimated 
from the diffusion constant assuming that the length in 
which the skyrmions effectively interact with the nucleus 
is about one lattice constant a. In this way r ~ a 2 /2D 
and 



To 1 = 
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Nf s (ln n) 



(3) 



Now, in order to make the temperature dependence 
of expression Q explicit, we need to know how the 
skyrmion diffusion constant changes with the tempera- 
ture. From the semiclassical point of view the study 
of the skyrmion dynamics was performedi^i^ starting 
from a generalized non-linear sigma model, which has the 
skyrmion as the lowest lying charged excitation 2 -. Within 
this approach the origin of the non trivial skyrmion dy- 
namics can be understood in terms of its interaction with 



the thermal bath of spin waves. This kind of interaction 
can be correctly treated using the well known collective 
coordinate method^. This technique promotes the cen- 
ter of mass of the skyrmion to a true dynamical variable 
and, as a final result, provide us with a Hamiltonian in 
which the skyrmion momentum is coupled to the gen- 
eralized spin wave momentum. The effective equation 
of motion for a single skyrmion was obtained using the 
Feynman- Vernon functional approach^ (by tracing out 
the magnons' degrees of freedom) and corresponds to a 
Brownian particle with temperature dependent damping 
and diffusion constants. The explicit temperature depen- 
dence of the skyrmion diffusion coefficients^ is given by 



D = DT 3 c X p(-2g/T), 



(4) 



where the Zeeman gap (g = g*/iBB) will be measured in 
Kelvin. 

In expression (JSJ not only the diffusion constant but 
also the number of free, or moving, skyrmions are tem- 
perature dependent quantities. In fact, as it was pointed 
out by Nederveen and Nazaro-^i 7 -, the donors situated on 
a layer located at a distance d from the 2DEG (inside the 
potential barriers) generate a random attractive poten- 
tial (disorder) in which the skyrmions start to localize 
below some temperature. For higher temperatures the 
number of skyrmions will be only determined by the de- 
viation of the filling factor from v = 1. Therefore, around 
any additional electron or hole (whose number is No) in- 
troduced in the completely polarized state, a new spin 
texture will be formed. In this case the number of free, 
or moving, skyrmions (Nf s ) exactly coincides with iVo, 
which is related to the electronic density of the 2DEG no 
by No = n S\l - v\. 

As it can be seen our model assume that as the tem- 
perature is lowered the number of free skyrmions de- 
creases. At first sight if the unbinding of skyrmions from 
the disordered centers is a thermally activated process, 
the number of free spin textures will be roughly given by 
N s f = noS |1 — v\ exp (-U/T), where U is some average 
value of the attractive potential induced by the disorder. 
However, this simple assumption implies that all parti- 
cles localize at the same temperature, in contradiction 
with the different patterns observed in the Knight shift 
measurements for different values of v. In order to be 
more realistic in modeling the disorder, a certain degree 
of correlation between them should be included. The 
main effect of this correlation is nothing but a distribu- 
tion of potentials with different depths in such a way that 
as the temperature decreases the skyrmions start to lo- 
calize in the deepest wells and this process continues until 
the "last skyrmion" localizes in the shallowest well. If this 
distribution is assumed to be Gaussian around some spe- 
cific value U , the number of potentials with depth U can 
be written as 

n(U)=Afexp(-(U-U ) 2 /A 2 ), (5) 

where the constant J\f is related to the total number of 
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donors localized in the barriers between the GaAs quan- 
tum wells. As all the electrons in the 2DEG arise from 
the donors their number should be equal, i.e., 



n Q S = N 



dUe 



-{U-U ) 2 /l 



= ^A(l + er/(Uo/A)), (6) 

where erf{x) is the error integral^. 

Now, using expressions £0 and © the temperature 
dependence of the transversal relaxation time (pj can be 
generalized to be a mean value over the potential distri- 
bution as 



dUe -{U-U f/^ e -U/T_ (?) 



Here, the value of U m in is related to the total number of 
skyrmions in the 2DEG through 



no\l-v\S = N 



dUe -(U-U f/A 

1 + erf ^ Umin ~ U ° 



(8) 



Notice that as the number of skyrmions increases 
with 1 1 — v\, the value of U m %n should decrease in order to 
accommodate all skyrmions in the localization centers. 

Finally, integrating Eq. Q and using the expressions 
© and O, it is possible to show that 
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where (j> is given by 
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1 + erf 
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2T 



Uo 
A 



(9) 



and erf~ 1 (x) is the inverse function of erf(x). 

As it can be checked expression © is a non mono- 
tonic function of the temperature. In fact in the high 
temperature limit, where the localization effects are neg- 
ligible, the NMR linewidth (T 2 _1 ) goes to zero in agree- 
ment with the motional narrowed effect observed in the 
experiments. If we look in the opposite direction (low 
temperatures) where the slow motion of the skyrmions 
tends to increase the linewidth, the exponentially small 
number of free skyrmions cancels out this effect and the 
linewidth decreases again. This is precisely the frozen 
limit found in the NMR measurements. Therefore the 
interplay between skyrmion dynamics and localization 
induced by disorder in the system reproduces some of 
the major features of the transversal relaxation time near 
v = 1, It is worth pointing out that depending on the 
disorder strength the 2DEG, away from v = 1, can evolve 
to different phases^ which in principle can influence the 
linewidth measurements. However, these physics was not 
included in our very simple model. 



III. SKYRMION IN THE PRESENCE OF 
DISORDER 

At this point nothing is left but to estimate the value 
of U$ due to the impurities localized outside the quantum 
well. This can be done computing the change in energy 
of an isolated skyrmion due to the presence of a posi- 
tive charged ion. In order to do that the skyrmions will 
be described by an effective generalized nonlinear sigma 
model^ in terms of an unit vector field m(r) associated 
to the electronic spin orientation. The Lagrangian den- 
sity which describes these objects in the presence of an 
external magnetic field B can be written as 



where 



Co = T(m) - V(m) 



T(m) = ^ A(m) ■ d t m, 



(10) 



(11) 



and 



V(m) 



p s (Vm) 2 + g*pp B m ■ B 
1 ' .n l(r) q(r') 



dr'- 



r — r 



(12) 



In the kinetic term iflOl) p denotes the electronic density 
and A[m] corresponds to the vector potential of a unit 
monopole, i.e., djA k = rrii. On the other hand in the 
potential energy density \V2l . p s is the spin stiffness, g* 
is the effective Lande factor, p = l/(2nl 2 ) is the uniform 
electronic background density and e is the dielectric con- 
stant of the background semiconductor. The deviation 
of the physical density from the uniform value p deter- 
mines the skyrmion density q(r), which can be explicitly 
written as 



q(r) = —e M m(r) ■ (<9 M m(r) X d v m{r)) 7 (13) 

and whose spatial integral is the topological charge. 

Instead of directly trying to solve l(Tol) - l|T2")l . Sondhi 
et al& considered the soliton solution of the nonlinear 
a model obtained by Belavin and Polyakov^ with a 
fixed size A, which is set by the competition between 
the Coulomb and Zeeman interactions included in {TJJ. 
Therefore, within this approach, the static solution of 
(IK)fl can be written as 



rr, 
m x{y) 



4Ar cos 9 (sin ( 
r 2 + 4A 2 



4A 2 



4A 2 ^ 



(14) 



which describes a skyrmion of unit topological charge 
localized at the origin. The skyrmion size is given by 
A = 0.558^(5 Ihigl) -1 / 3 , where lo stands for the mag- 
netic length and g = gp, B Bje. In this approximation the 
skyrmion density is given by q{r) = 4A 2 /7r(4A 2 + r 2 ) 2 . 
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r <MP 



Figure 1: Effective potential in if as a function of ro ( the 
skyrmion-donors distance in the the 2DEG). The continuous 
line corresponds to d = 900A and the dashed line to d — 
1800 A. 



As pointed out earlier, the disorder effect in the 2DEG 
is related to the donors situated in a layer whose dis- 
tance from the 2DEG is d. In order to include those ef- 
fects in the skyrmion dynamics, we add an extra term to 
the the Lagrangian density 11 Oil . If we consider that the 
skyrmion density q(r) is coupled to the disorder poten- 
tial via Coulomb interaction, our model for the skyrmion- 
impurity system is given by the Lagrangian density 



C = L a 



2 q(r) 



e|d 2 + (r-r ) 2 | 



(15) 



where e is the electron charge and ro denotes the donor 
coordinate at the 2DEG plane. Once we are assuming a 
very simple model to describe the system, we will not con- 
sider the fact that the disorder potential can be screened 
by the 2DEG as pointed out by Efros et al—. 

>From expressions l(Tflll - lfT2|l and l(T5|l . the energy func- 
tional for a single static skyrmion can be written as 



E 



V[m] 



2 q(r) 



e|d 2 + (r-r )2| 



dr 2 



(16) 



where V(m) is the interacting potential related to the 
Lagrangian density dlOfl - J12f) and the integral of the sec- 
ond term is performed over the whole plane containing 
the 2DEG. If we assume that the presence of the dis- 
order does not modify the skyrmion form and size, the 
energy of the system will be given by Equation 11 fill with 
m(r) = m°(r). Therefore, the change in energy of a sin- 
gle skyrmion due to the disorder potential is given by 
U eff (r ) =E-E , where E = J 7(m°(r, X))dr 2 . 

The specific form of U e ff as a function of rg (the 
skyrmion-donor distance in the plane of the 2DEG) is 
illustrated in Fig. Q]for A = 1.2Z (the skyrmion size at 
B = IT , see appendix). As it can be seen, the minimum 
in energy corresponds to the case where the impurity 



position in the 2DEG plane exactly coincides with the 
skyrmion center. Therefore, the parameter C/ in Eq. © 
can be estimated as the value of U e ff at r = 0. If we 
use a typical value of the distance from the donors to the 
2DEG as d ~ 1.800 Alft, we conclude that U Q ~7K. It 
is interesting to notice that the value of the estimated Uo 
is of the same order of the disordered localizing potential 
reported in Ref . for the case in which v w 1 . 



IV. RESULTS AND DISCUSSION 

At this point we can turn or attention to the linewidth 
behavior as a function of the temperature. Fig. shows 
a plot of equation © as a function of temperature for 
different values of |1 — v\ with t/ = 7K, A = 2 K and 
g = 2 K (the Zeeman gap at B = 7T), One can see that, 
T 2 _1 has the non-monotonic behavior observed in the OP- 
NMR experiments^^ going from the motional-narrowed 
regime at relative high temperatures to the frozen limit 
as the temperature decreases. The agreement between 
the theory and the experimental data reported in Ref. 
S suggests that the relevant gapless mode that induce 
short nuclear relaxation times Ti near v w 1 is precisely 
the translational mod o 1 ^' 14 . On the other hand, it can 
be seen that as | 1 — v | increases the maximum value 
of 1/T 2 also increase. On physical grounds this can be 
expected, because as the number of skyrmions grows, the 
possibility of finding free skyrmions which effectively in- 
duce nuclei dephasing, increases. Although for the 
measurements^ this kind of behavior is almost unclear, 
in the v w 1/3 situation it can be neatedly observed^. 

For Uq equal to the potential depth estimated from the 
isolated impurity case, the peak position -indicating the 
change in the dynamical regime- is around 1.7 K which is 
very close to the reported experimental data£. It is possi- 
ble to show that if Uq changes from 4 K to 8 K, the peak 
position varies from 1 K to 2.2 K remaining within the 




Figure 2: T 2 1 (FWHM) as a function of temperature for 
Uo = IK, A = 2K, g = 2K and different filling factor 
values. 
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reported results^. As it can be seen in Fig. the peak 
position is almost constant as the filling factor changes, 
which is in contradiction with the experimental results for 
v ~ 1» but in agreement with the v s» 1/3 caseiS. This 
discrepancy could be possibly associated with the non in- 
teracting skyrmions approach we have employed, that is 
more appropriate for the small skyrmion near v = 1/3. 
This conjecture can be tested by applying an in-plane 
magnetic field, which does not change the filling factor 
but increase the Zeeman energy making the skyrmion 
smaller. 

It is worth pointing out that our model for the 
linewidth of the NMR signal is in agreement with the 
scenario suggested for the observed tilted plateau in the 
in the Knight shift measurements for the 2DEG when 
v w 1. In this case, the existence of the tilted plateau 
is also related to the localization of the quasiparticles in- 
troduced in the 2DEG as the filling factor deviates from 
v = 

Summarizing, we have developed a very simple theory 
which describes the evolution of the NMR line profile of 
the 71 Ga nuclei coupled to the 2DEG (near v = 1) as a 
function of the whole available temperature range which 
seems to be in fairly good agreement with the experimen- 
tal data. 
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VI. APPENDIX 

In the table below, we show the energy scales for the 
quantum Hall system at v = 1 in Kelvin. In all expres- 
sions, the magnetic field is measure in Tesla. For GaAs 
quantum wells e w 13. 



Energy scales 


(K) 


e c e 2 /d 
9 9*PbB 
Ps e c /(16V2^) 


50.4O/B 
0.33B 
1.25VB 



V. ACKNOWLEDGMENTS 

AVF would like to thank A.R. Pereira and M.O. Goer- 
big for helpful discussions. The authors kindly acknowl- 



Therefore the skyrmion size A = 0.558Zo(g |hi<?|) 1 / 3 
calculated as described in the text can be written as 

A = 3.08 WB In (5.95 x IQ-^^fBj ) , where l = 

^Hc/eB = 256/VsA. 
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